Virial mass in DGP brane cosmology 



H. R. Sepangi* and S. Shahidi^ 

Department of Physics, Shahid Beheshti University, Evin, Tehran 19839 Iran 

August 19, 2009 



Abstract 

We study the virial mass discrepancy in the context of a DPG brane-world scenario and show 
that such a framework can offer viable explanations to account for the mass discrepancy problem. 
This is done by defining a geometrical mass M that we prove to be proportional to the virial 
mass. Estimating M using observational data, we show that it behaves linearly with r and has a 
value of the order of M200, pointing to a possible resolution of the virial mass discrepancy. We 
also obtain the radial velocity dispersion of galaxy clusters and show that it is compatible with 
the radial velocity dispersion profile of such clusters. This velocity dispersion profile can be used 
to differentiate various models predicting the virial mass. 

1 Introduction 

The past decade has been marked by the advent of an avant garde school of thought which has tried 
to address a number of important problems and observations in standard General Relativity and 
cosmology, namely the hierarchy problem; the great disparity between the fundamental forces of 
nature, the observation that the universe is undergoing and accelerated expansion relating to dark 
energy and the problem of galaxy rotation curves relating to dark matter. This school of thought 
is based on the assumption that our 4-dimensional observable universe, the brane, is embedded 
in a higher dimensional space, the bulk, which has the geometry of an AdS space and to which 
gravitons can scape but ordinary matter cannot. The AdS nature of the bulk space would cause 
gravity to become localized around the brane and would modify the gravitational potential at short 
distances. Such a scenario pQ, proposed by Randall and Sundrum (RS), has been able to account 
for the hierarchy of the fundamental forces with great success and has been generating a myriad of 
other scenarios and variations. Shortly after, in a seminal work [2], it was shown how to project the 
Einstein field equations, assumed to hold in the bulk, onto the brane. An unprecedented number of 
works utilizing this idea have been appearing ever since O 0]. 

In an effort to relax the restriction of an Ads bulk and hence allowing gravity to penetrate large 
distances, Dvali, Gabadadze and Poratti (DGP) [5j [6] proposed an alternative model in which the 
influence of gravity on the brane is accounted for by including an induced 4-dimensional Einstein- 
Hilbert term to the full action. In this model, in contrast to the standard RS model, gravity is 
modified at large distances. For a comprehensive review of this model see [TJ. The cosmological 
implications of this model was investigated in [H] where it was shown that the Friedmann equation 
on the brane has two branches, both reducing to the usual FRW equation at the small Hubble radius 
limit. However, the important discovery was that one of these branches predict a self-accelerating 
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universe at late times, consistent with the observation that our universe is undergoing an accelerated 
expanding phase. So much for the success of the DGP model, a word of caution is in order; the 
theory predicts the existence of ghost-like excitations. Many scenarios have been undertaken to 
explain away such ghosts, but as yet no satisfactory solution exists. The interested reader should 
consult [HUH] for further insight. In this paper we do not discuss such excitations since our aim lies 
in studying the virial mass discrepancy in DGP models. 

One of the interesting problems in cosmology is the calculation of the mass of cluster of galaxies. 
In recent years, our ability in performing precision measurements in observational cosmology is 
developed to such an extent that we can obtain accurate values for the mass of individual galaxies 
and their velocities. We may therefore find the total mass of cluster of galaxies in two ways, adding 
up the masses of individual galaxies, or do it statistically and use the virial theorem. Since these 
methods represent two aspects of the same thing we must obtain the same result. However, almost 
in all clusters the virial mass is 20-30 time greater than M, the mass obtained by adding up the 
individual masses. Also by Newton second law we know that the mass of a galaxy is proportional to 
rv tg where vt g is the tangential velocity of a test particle located at the distance r from the center 
of the galaxy. As observations have shown, the tangential velocity of a test particle remains nearly 
constant at large distances from the center of the galaxy. This is, of course, in contradiction to 
what Newtonian gravity predicts. One way around this is to postulate the existence of dark matter. 
However, there are geometrical approaches to address this problem, namely to use modified Einstein 
field equations, as is done in brane- world models \\.\\ I12j or in modified gravity |13t I14j. 

In this paper, we use the DGP model discussed above to explain the the virial mass discrepancy 
in a geometric manner. As we shall see later, the linearly increasing behavior of the virial mass with 
distance can be explained by taking into account the extra terms that appear in the field equations 
which in turn originate from the bulk geometry. In order to do this we must have a procedure to 
obtain the total mass of the cluster. We use Jean's equation and the observational data to obtain 
the total mass distribution of clusters. We also obtain the radial velocity dispersion in this model 
which suggests an alternative way to obtain the virial mass of clusters and can be used to explain 
the observational data. The question of the flat rotation curves of individual galaxies requires a 
separate undertaking and will be dealt with in a future work. 



2 Einstein equations on the brane 

Let us start with the standard DGP action [6] 

/ d 5 x^iJZ + ^ [ d 4 Xy /=q~(R - 2A) - m\ [ d 4 x^K + S m (q^) + S B (g AB 
JM 1 JdM JdM 

(1) 



s = ^ 



where gABi 1Z and Sb are the metric, Ricci scalar and matter action of the bulk and q^, R and 
S m are those of the brane with A being the brane cosmological constant, and m\ (m|) is the bulk 
(brane) Planck scale. The third term is the Hawking-Gibbons boundary term [15] and K is the 
extrinsic curvature. After varying the action and denoting the extra dimension by y, we obtain the 
Einstein field equations in the bulk 



8^8 B (Tpuv - m|A<^)<%) + f AB , (2) 
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where Tab {T^p) is the bulk (brane) energy momentum tensor for which a perfect fluid form is 
assumed 



T A B = diag(-p B ,p B ,p B ,p B ,p 5 ), 
T*i = diag(- Pb,Pb,Pb,Pb,fy, 

and for the bulk metric we take 

ds 2 = -e uir ' y) dt 2 + e^ r ' y) dr 2 + r 2 (d6 2 + sin 2 6 d^ 2 ) + e x{r ' y) dy 2 . 
The induced metric on the brane is simply 

ds 2 = - e uo ^ dt 2 + e^°W dr 2 + r 2 (d6 2 + sin 2 9 d^ 2 ) , 
where z^o(r) = ^faO) etc. Assuming that the brane is at y = and noting that [16] 



d 2 / 

dy 

df\ (dh 



- r L = f + 2f8(y), 



dy 1 \dy 



fh, 



(3) 
(4) 

(5) 

(6) 

(7) 
(8) 
(9) 



where / = and /, h are arbitrary functions of \y\, we obtain the following field equations in the 
bulk 

m\ 
4^2 
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where a prime represents derivative with respect to r. Since only 5(y) can contribute to the brane 
equations we obtain the Einstein field equations on the brane 

(15) 
(16) 
(17) 
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where we have defined the induced energy density U(r) and pressure V{r) as 

U{r)=mlfi\ y=0 e~ x \ (18) 
V{t) = -ml u\ V=0 e- X °. (19) 



3 The Virial theorem 

To obtain the virial theorem in the context of the model discussed above we use the tetrad formalism 
by defining the following frame of orthonormal vectors [T7] 

ef=e^5% eW=e^J, ef = r8% ef=r^6 5% (20) 

where q^eff e^ = rf a W) and the tetrad indices are surrounded by parenthesis. The 4- velocity 
of a typical galaxy with v^v^ = — 1 is written as 

»W=AW a = 0,1, 2, 3. (21) 

Let us start with the Boltzmann equation in tetrad formalism. If f(x^, v^) represents the distribu- 
tion function of galaxies, supposed to be made of identical and collisionless point particles, we have 

03 EES] 

v (a) P ^/ ■ (a) (b) (c) d f _ Q (22) 

where 7m( c ) = e p a ]^ e (b) e (c) are ^ ne Ri cc i rotation coefficients. Using the brane metric © and 
assuming that f(x fl ,v^) depends on r only, the Boltzmann equation becomes 

df_ 

■ 0. (23) 

/ \ U V(p U U0 J 

where we have defined 

«(°) = Vt , V® = V r , V (2) = Vg, W (3) = V„. (24) 
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Since our metric is spherically symmetric, the coefficient of cot must be zero in (]23p . Multiplying 
equation (|23|) by mv r dv, where dv = ^dvrdvgdv^ is the invariant volume element in the velocity 
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space and m is the mass of the galaxy, and integrating over the velocity space and assuming that / 
vanishes sufficiently rapidly as the velocities tend to ±00, we obtain 

^ [P ( V r)} + \P K V t) + ( V r)} T < ~ P i( v2 e) + ( V D ~ 2 ( V r)} = °> ( 25 ) 

where p is the mass density and ( ) represents the average values. Multiplication of equation (|25|) 
by 47rr 2 and integration over the cluster of galaxies yield 

\p [{vl) + {v 2 ) + (vl)] r 2 dr + \ j\nr 3 p [<«?) + <, r 2 >] ^dr = 0. (26) 

We can also write equation (|26|) in the form 

2K = \j\^p[(vl) + (vl)] d -^dr, (27) 

since the total kinetic energy of galaxies is defined as 

K= [ R 2tt P [{v 2 r ) + {v 2 ) + (4)] r 2 dr. (28) 
J 

To obtain the virial theorem in our model we must express the energy-momentum tensor components 
in the terms of the distribution function. This is done according to 

T - 
- L fj,i> - 



^ fmv fJ/ v u dv, (29) 
which leads to 

Pb = p{vt), Pb = p(vr) =p(«|) =p(«J). (30) 
Adding equations (fTST) . (fl6j) and twice of (fTTl) yields 



mfe^ (7 - ^ + Y + f ) = ^ ^ + \ ~ U{T)] ~ miA ' (31) 

where we have defined {v 2 ) = {v 2 ) + {v 2 ) + {vq) + {v 2 ^). For the cluster of galaxies we may assume 
that p(r) and u(r) are small so that the quadratic terms in equation (|3ip do not contribute. Also, 
the velocity of galaxies are much smaller than the speed of light, so we can set {v 2 ) , {vg) , {v 2 ) <C 
{v 2 ) w 1 [TT]. Taking these approximations into consideration, equation (J3T1) is reduced to 



P = {r 2 u' ) + 2m|A - [3P(r) - W(r)] . (32) 

Multiplying equation (|32|) by r 2 and integrating from to r yields 

12 1 

m?r 2 i/„ M (r) + -mi Ar 3 J\f(r) = 0, (33) 

where 

/•r 

M(r) = 4vr / p/dr', (34) 

JO 
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and 



rr 

J\f(r) = 4vr / [3V(r') - U{r')] r' 2 dr'. (35) 
Jo 

Again, multiplying equation ()33[) by dM ^ an d integrating from to R, we finally obtain the 
generalized virial theorem in a DGP senario 

W + 2K + \kI + W B = 0, (36) 



where 



1 f R M(r 



W = -—^ / — ^dM(r), (37) 



f 

'3 JO 
and 



W B = --^ / prjV(r)dr, (38) 
2m, 



/= / r 2 dM(r), (39) 
JO 

is the moment of inertia of the system. Without the last term, this would constitute the usual virial 
theorem with a cosmological constant, first derived by Jackson [T7] using the Boltzmann equation 
(|22[) into which the metric of the space-time is substituted. In order to obtain a relation between 
the virial mass and the extra term M(r) which has its origins in the bulk, we define the following 
radii [H] 

M 2 

R v = ttj \ , (40) 

y jRMirl dM ( r y 

n f R r 2 dM(r) 

- S?T' < 41 > 

K = -8^fC' (42> 

where R v is the virial radius and 1Z is the radius defined by the extra term J\f. By defining the 
virial mass as 

1 M 2 

2K = ^2^^ (43) 



and using the relations 



87rm3 R v ' 



1 M 2 

^"wp/ I = MR - < 44 > 



the generalized virial theorem (|36p can be written as 



My\ 2 = SnrnjA R 2 R V / Af\ 2 f R, 
Ml 3 M \M \ n ' ' 
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The contribution of A to the mass of the galaxy is several order of magnitude smaller than the 
observed mass. Also M v is much larger than M for most galaxies. Therefore, we can neglect the 
unity and the term involving the cosmological constant in equation (|45H . The virial mass in our 
model is then given by 



M v {r)~N(r)^. (46) 

As can be seen, the virial mass is proportional to an extra term stemming from the global bulk 
effects. 



4 Estimating J\f(r) 

In order to estimate J\f(r) we must solve the Einstein equations for V(r) and U(r) and obtain J\f{r) 
from (I35p . However there is a simpler way of doing this which we will follow. First, consider the 
conservation of the right-hand side of Einstein equations (|15p -(|17 p 

W 



p' b + [v'+f i 

V '^~ 2 WT^TWTuy (47) 

U = 0. (48) 

This means that we only need to calculate V(r). In most clusters the majority of the baryonic mass 
is in the form of intra-cluster gas. Taking this assumption into consideration and using equations 
(|34"|) and (f35|) . we obtain an expression for the total mass of the cluster 



dM tot 



4irp g r 2 + 12-KVr 2 . (49) 



dr 

Another expression can be obtained from Jean's equation 

|W]+P,(r)§ = 0, (50) 

where <3?(r) is the gravitational potential. We also assume that the gas is isotropically distributed 
inside the cluster so that the mass- weighted velocity dispersion in the radial and tangential directions 
are equal; a r = ag^. Assuming that the gravitational field is weak so that 3>(r) satisfies the poisson 
equation 2m|V 2 <I > « pt t, the Jean's equation is reduced to 

dp g {r) 1 M tot 



dr 87rmi 



: Pg(r), (51) 



where we have used the relation p g = p g a 2 . As has been shown in [TO], the gas density p g can be 
fitted to the observational data by the following radial distribution 



■2 



p g {r)=p [\ + 71 J , (52) 

where r c is the core radius and /?o an d P are cluster independent constants. For most clusters /3 > | 
[19] and therefore, in the limit r r c considered here, the gas density distribution can be written 

as 

\ -3/3 o 
r \ 2 



P 9 (r)=Po(-) , /5>3- (53) 
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Moreover, we choose the following equation of state for the intra-cluster gas |19j 



(r) = KIl (r) (54) 

where [i = 0.61 is the mean atomic weight of the particles in the cluster gas and m p is the mass of 
proton. With these assumptions, equation (|5T|) reduces to 

3£- T 

M tot (r) = 8<iTm 2 ^^f3r. (55) 
fxm p 

The contribution of the gas density to the total mass of the cluster is very small at the boundary of 
the cluster, i.e. where r 3> r c , as can be seen from equations (I53D and (I49p . We may now calculate 
V{r) from equations (|49p and (|55p . taking the approximation above 



> k D ^ 1 



V(r) = 2m 2 3 ^^(3 -j, (56) 



and using equation ([35]) to obtain 



3fr T 

M(r) = M tot {r) = %irnil—^l3 r. (57) 

This equation is obviously an approximation since the contribution of the baryonic mass is neglected. 
However this is no cause for concern since the value of the baryonic mass of the clusters is about 
3 orders of magnitude smaller than its total mass. As can be seen, M{r) oc r and since the virial 
mass is proportional to M and the latter is proportional to r, this could offer a possible resolution 
to the virial mass discrepancy in the context of DGP brane worlds. 

To estimate the value of Af, we first note that 87rm 2 = GT^ = f G" 1 where G N is the gravitational 
constant [Sj. A typical value of the temperature of a cluster gas is k B T g ~ 5 keV [19] . The virial 
radius of the cluster of galaxies is usually assumed to be r2ocb indicating the radius for which the 
mass density of the cluster is about P200 = 200/0 C , where p c = 4.6975 x 10" 27 h^ kg/m 3 . The virial 
mass of the cluster is then estimated as My = M200 = M(r < r2oo)- We can therefore define the 
maximum extension of the M{r) mass to be the radius at which V = /9200- From equation (|56p we 
have 

r max = 4.28/3^ h£ (^J^) 2 Mpc. (58) 



k D T n \i 



Finally, we can estimate M(r) from (I57p 

J\f(r) = 32.72 x 10 14 /35 ( "f/ g TA ) M . (59) 

\5 KeV J 

This is in agreement with observational values for the virial mass of clusters [19J. 



5 Radial velocity dispersion 

Another important observational quantity is the radial velocity dispersion which plays an important 
role in estimating the virial mass of the clusters. As is well known, the simple form a 2 = B j(r + b) 
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Fi gure 1: Radial velocity dispersion for the NGC5813 cluster with f3 — 0.766. For this cluster r2oo — 0.87Mpc, 
r c = 25Kpc, k B T a = 0.52.KW and we have set Ci = 4.2, 2.9, 1.4 x 1O" 8 M and C 2 = 0.0390, 0.0260, 0.0130 M^/Kpc 4 
for the top, middle and bottom curves respectively. 



can be used to fit the observational data [20J. The DGP model can provide an expression for the 
radial velocity dispersion which we will derive in this section. 

Adopting the approximations used after equation (j3Tj) . we may write equation ([25]) as 



-|-;(^r) + 2^o = 5 



(60) 



where we have assumed that the velocity distribution in the cluster is isotropic, so that (i> 2 ) 
(vfy + (vq) + (yfy = 3 (y^ = 3of . Moreover, from the Einstein field equations we have 



ml 



=3V(r)+p(r). 



Integrating, we obtain 



2 2 / 
m 3 r u 



(61) 



(62) 



where C\ is some constant. The differential form of the radial velocity dispersion can be obtained 
from equations (f60j) and (f62|) 



M{r) M(r) Ci 



Now, using expressions (f53"|) and (|57p for p(r) and J\f(r) and by virtue of equation 



a, 



2 _ ^B^9 



fim p 12(/3- l)(3/3- l)m 



3 V r c 



-3/3 
- I + 



1 

- + 



Co 



2(3/3 + l)mg r 2/711/90 



3fl 



(63) 
we obtain 

P^l. (64) 
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For (3 = 1 we have 



a 2_KT g | Porglnr | C 3 | C 2 ( r\ 3 
finip 8m| r r 2m|po \ r c/ 

Our expression for can therefore be used to fit the observational data. In figure [1] we have 
plotted the radial velocity dispersion for the cluster NGC5813. This cluster has (3 = 0.766 and 
ksTg = 0.52KeV |19j . and the radial velocity dispersion is about 240/cm/s |21j . We see that the 
radial velocity dispersion (j64"|) is compatible with the observed profiles [211 [20]. Since an expression 
for the radial velocity dispersion can be obtained by other theoretical methods which explain cluster 
discrepancies, such a relation can be used to differentiate them 122]. 



6 Discussion 



In this paper we have considered the virial mass in the frame work of DGP brane worlds. The 
resulting field equations on the brane have an additional term which is due to the geometry of the 
extra dimension and can be associated with a geometrical mass. The virial theorem was obtained by 
the use of Boltzmann equation, assuming that galaxies in the cluster are point-like, non-interacting 
particles. We showed that the resulting virial theorem has an additional potential term due to the 
extra dimension. The virial theorem has also been exploited in other brane- world models [111 [T2"] to 
explain the virial mass discrepancy. The advantage of the DGP model however is in the explanation 
of the self accelerating phase of the universe without resorting to dark energy in a consistent manner. 
We obtained the virial mass of clusters from the virial theorem and showed that it is proportional 
to the geometrical mass of the model. The behavior of the virial mass was investigated in section 4 
and showed to be a linear function of the distance. To estimate the geometrical mass M we needed 
another relation in addition to the Einstein field equations to close the system of equations. Such a 
relation was obtained by using the Jean's equation for clusters, assuming that clusters have spherical 
symmetry and are in thermodynamic equilibrium. 

The solution presented in this work offers a possible explanation to the question of the virial mass 
discrepancy. The radial velocity dispersion profile of clusters was also obtained, having two arbitrary 
constants, one with a coefficient decreasing with r and the other increasing with r. This allowed us 
to use this expression for any cluster with (3 > |. We used the observed radial velocity dispersion for 
the cluster NGC5813 as an example to show that our model can account for the velocity dispersion 
of clusters. In addition to explaining the observational data, the velocity dispersion profile can be 
used to study the various aspects of models predicting the virial mass. 
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